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Abstract 

' Data on pp annihilation in flight into Tv^^Ti^'ri are presented for nine beam momenta 600 to 1940 MeV/c. 

^ The strongest four intermediate states are found to be /2(1270)r;, a2(1320)7r, arj and ao(980)7r. Partial wave 

Jj~^. analysis is performed mainly to look for resonances formed by pp and decaying into tt^tt"?; through these 

intermediate states. There is evidence for the following s-channel 7 = resonances : two 4"*"^ resonances 
with mass and width (M, T) at (2044, 208) MeV and (2320 ± 30, 220 ± 30) MeV; three 2++ resonances at 
. (2020 ± 50, 200 ± 70) MeV, (2240 ± 40, 170 ± 50) MeV and (2370 ± 50, 320 ± 50) MeV; two 3++ resonances 

at (2000 ±40, 250 ±40) MeV and (2280 ±30, 210 ±30) MeV; a 1++ resonance at (2340 ±40, 340 ±40) MeV; 
j^'! and two 2"+ resonances at (2040 ± 40, 190 ± 40) MeV and (2300 ± 40, 270 ± 40) MeV. 
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1 Introduction 

> , 

, The Crystal Barrel detector is being used to make a systematic study of the mass region 1960 to 2410 MeV in 

OP ' pp annihilation in flight at LEAR, with p beams of momenta 600 to 1940 MeV/c. The objective is to study 

, resonances in the formation process, i.e. the s-channel. Here we study data in pp — ^ ir'^'K^rj for resonances 

decaying to a2(1320)7r°, /2(1270)77, /o(1500)77, /o(980)77, ao(980)7r° and a-q. We use a to denote the broad tttt 
S-wave amplitude up to 1860 MeV. The present results have been presented briefly in the form of a letter 
[1] and here we give full details of the experimental techniques and analysis. Further studies of 777777° have been 
presented elsewhere [2,3], and work is in progress on other channels such as 37r" and 7r"7r"77'. 

From earlier work, it is known that the mass range we explore contains many resonances [4] ; a detailed study 
of pp — TT^TT"*" using a polarised target has provided much of the current evidence [5,6]. The /4(2050) is well 
known, and from the quark model of meson resonances one expects that it will be accompanied by /s and /2 
resonances close-by in mass. We shall indeed provide evidence for these resonances and a further one with 
quantum numbers J^^ = 2 and similar mass. At higher masses, towards the top of the LEAR range, there 
has been evidence for /4(2300) and /2(2340) [4], and it is anticipated from the Veneziano model [7] that there 
is likely to be a tower of resonances around this mass. We shall provide evidence for such states with quantum 
numbers 4+, 3+, 2+, 1+ and 2". 

These resonances are anticipated qq states. This mass range is also likely to contain glueballs with quantum 
numbers and 2++, predicted in the mass range 2000-2400 MeV by various theoretical models [8,9,10]. 
Hybrids may also be present. Decays of these exotic resonances to 77 and a seem to be favoured in /o(1500) 
decay [11], charmonium decay and J/^P radiative decays [12]. Hence the 77/2(1270) and t/ct channels are of 
particular interest. 

The layout of this paper will be as follows. In section 2, the procedure for data processing and event selection 
is outlined; the data are presented and their gross features are discussed. Section 3 gives the formalism used 
for the partial wave analysis. Section 4 gives the results for partial wave amplitudes. Then, in Section 5 we fit 
partial waves to resonances. Finally, Section 6 provides a summary. 
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2 Experiment and Data Processing 



The data were taken at LEAR by the Crystal Barrel Collaboration, using a trigger on neutral final states at nine 
beam momenta from 600 to 1940 McV/c. An average of 9 x 10^ triggers were taken at each momentum. The 
detector has been described fully in an earlier publication [13]. 

A liquid hydrogen target 4.4 cm long is surrounded at increasing radii by a silicon vertex detector, a multiwire 
chamber for triggering, a jet drift chamber to detect charged particles and finally 1380 Csl crystals to detect 
photons. The present data were taken with a trigger demanding a neutral final state. For this purpose, the 
silicon vertex detector, multiwire chamber and jet drift chamber were used simply to veto charged particles. 

The barrel of Csl crystals covers 98% of 4tt solid angle. Crystals are 16 radiation lengths long and point towards 
the target. The angular resolution is ~ ±20 mrad in both polar angle and azimuth. The detection efficiency is 
high for photons down to energies below 20 MeV. The energy resolution AE is given by AE/E = 0.025/E^/^, 
where E is in GeV. 

The incident p beam was pure and monoencrgctic with momentum spread Ap/p < 0.1%. Incident antiprotons 
were defined by a coincidence between a small proportional counter P and a 5mm diameter silicon counter, Si. 
Two veto counters, 20 cm downstream of the hydrogen target, were used to provide a first level trigger P.Si.V 
identifying interactions in the target. The beam intensity was typically 2 x 10^ p/s and at times was twice this. 
The interaction rate in the target (excluding pp elastic scattering, where the forward p generally counted in the 
veto counter) was typically 3KHz. Of this, ^ 1 — 2% consisted of neutral final states, so the trigger rate for 
all- neutral events was 20-60 Hz. In order to filter out events which obviously fail to conserve energy, the total 
energy in the Csl crystals was summed on-line [14] ; a fast trigger rejected those events with total energy falling 
~ 200 MeV or more below that of pp annihilation. 

The absolute normalisation is derived from beam counts P. Si, target length and density, the number of 
detected events and a Monte Carlo simulation of reconstruction efficiency in the Csl barrel. Details of this 
normalisation are given in a paper on the tt^tt" final state [15]. A dependence of the reconstruction efficiency on 
beam rate is observed, and the normalisation has to be obtained from an extrapolation to zero beam rate. The 
normalisation uncertainty is estimated as ±3% from 1800 to 1050 MeV/c and increases to ±6% at 900 and 600 
MeV/c. Data at 1940 MeV/c were taken in separate, earlier runs, and have an estimated uncertainty of ±10% 
in normalisation. There is in addition an overall normalisation uncertainty of ±2.4% from the target length, 
common to all momenta. 

2.1 Data Selection 

A large number of alternative prescriptions have been examined for selecting events. At high momenta, one 
of the problems is that photons from decay sometimes merge into a single shower. Conversely, one shower 
sometimes splits into a primary shower and a nearby secondary shower, caused by Conipton scattering. The 
probability that this occurs is ~ 10% per photon. In early studies, an attempt was made to salvage ijn^n^ events 
from 57 or 77 final states. However, the gain in statistics was small 15%) and the penalty was an increase in 
backgrounds. Eventually, it was decided to retain only events containing exactly 6 photon showers. 



Beam Momentum 


Number of Events 


Reconstruction 


Cross Section 


(MeV/c) 




Efficiency (%) 




600 


20385 


26.3 


71.9 ±3.6 


900 


112476 


25.4 


83.2 ±4.9 


1050 


86238 


24.9 


78.9 ±2.3 


1200 


124581 


24.2 


68.6 ±3.0 


1350 


81454 


23.4 


54.4 ±2.3 


1525 


57714 


22.7 


56.5 ± 1.8 


1642 


65984 


21.9 


53.2 ±2.5 


1800 


71738 


20.8 


43.8 ±1.5 


1940 


75325 


19.8 


37.0 ± 3.7 



Table 1: Numbers of selected events, reconstruction efficiency and cross sections for pp — )• tt^^tt^t] with rj — > 77. 
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Data arc fitted kincmatically to a large number of physics channels: 43 for (4 — 8)7. In order to assess 
branching ratios to every channel and cross-talk between them, we generate at least 20,000 Monte Carlo events 
for every one of the 43 fitted channels, using GEANT. In the first approximation, events fitting the correct 
channel determine the reconstruction efficiency in each channel. Events fitting the wrong channel estimate 
the probability of cross-talk Xij between channels i and j. More exactly, we solve a set of 43 x 43 simultaneous 
equations containing on the left-hand side the observed number of fitted data events Di, and on the right-hand 
side reconstruction efficiencies and true numbers of created events Ni in every channel and terms allowing for 
cross-talk Xij between channels: 

A = €iN,+Y^ x.jNj. (1) 
The solution is constrained so that the numbers of real events, Ni, in every channel are positive or zero. 
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Figure 1: Cross section for pp ir^iT^ri with r] 77. 




This procedure is carried out for a variety of confidence levels (1, 5, 10, 20%) and using a wide variety of 
selection procedures. A choice is then made, optimising the ratio of signal to background. We find that this 
ratio is not very sensitive to confidence level over the range 5-20% for TyTr^Tr^ events. 

Among 6 — 8 photon events, the four largest channels are 47r°, 37r", ir^n^r] and tt^tt^w with w — >■ 7r°7. The 
relative branching ratios for these channels are roughly 1.1 : 1 : 0.4 : 0.4 at 1800 MeV/c. To select the w'^TT^r] 
channel, we demand exactly 6 photons satisfying a 7C kinematic fit with confidence level > 10%; events fitting 
377° with confidence level > 0.01% are rejected, and also those few events fitting tt^tt^t]', 7r°7y?y, n^rjr]' and 3r] 
with confidence level larger than that for 7r° 77^77. The Monte Carlo simulation shows that the worst backgrounds 
arise from wTr^Tr", (oj — > 7r''7) when one photon is lost, and from Att'^ events when two photons are lost. Residual 
backgrounds from these two processes are 1.5% and 0.8% respectively at 1800 MeV/c. Including other small 
backgrounds, the total is 3.0 ±0.3% at 1800 MeV/c. For lower beam momenta, the background increases slightly. 
At 600 MeV, the total is 4.0 ± 0.4 with the worst backgrounds from u;7r°7r° (1.7%), 47r° (0.9%) and looj (0.9%). 
Table 1 summarises numbers of selected events, the reconstruction efficiency and cross sections. Statistics at 600 
MeV/c are lower because most data were taken without the threshold cut on total energy in the trigger. The 
cross sections for the vyTr^Tr" channel are also shown in Fig.l. There are clear enhancements at low mass and 
around 2200-2300 MeV. Note that for a constant amplitude the cross sections should decrease steadily as the 
energy increases, see equn (23) below. 

Fig. 2 shows the confidence level (CL) distribution for data of beam momentum at 1.2 GeV/c. The slight 
peak at high confidence level arises from events where all particles emerge close to the beam direction, with 
the result that the vertex is poorly defined. We apply no cut on the coordinate of the vertex along the beam 
direction, so as to avoid biasing the data selection. The rise at low confidence levels is followed accurately down 
to 10% by the Monte Carlo simulation; it arises from overlapping showers in the Csl detectors. 

In order to illustrate the cleanliness of the 77 signal, we have made an additonal fit to 7r° 77^77. Fig. 3 then shows 
the mass distribution of 77 pairs in the vicinity of the t] peak for CL(7r°7r°77) > 0.1 with C L{-k^ tt'^ rf) > 0.0001 
at beam momentum 1.2 GeV/c. The 77 peak is well centred at the correct mass, 547.5 MeV and the background 
under the ry signal is compatible with that expected from the Monte Carlo simulation. 
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Figure 2: Distribution of confidence level for pp -k^-k^t] events at beam momentum 1.2 GeV/c. 
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Figure 3: Mass distribution of 77 pairs for Ci(7r°7r°77) > 0.1 with CL{'K^-K^rf) > 0.0001 at beam momentum 
1.2 GeV/c. 
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Figure 4: Dalitz plots for pp — > w^iT^r] at incident beam momenta 0.6 — 1.94 GeV/c. 



2.2 Features of the Data 



Fig. 4 shows Dalitz plots at the nine available momenta and Figs. 5 and 6 projections on to rjir mass and tttt 
mass. The most prominent feature of the Dalitz plot consists of a diagonal band due to /2(1270)7r°. There 
are weaker horizontal and vertical bands due to ao(980)7r and O2(1320)7r. The /2(1270)r7 signal grows with 
respect to a2(1320)7r as the beam momentum rises; this is a natural consequence of the increasing phase space 
for /2(1270)?7, whose threshold is at 1820 MeV. Very weak peaks are visible in the nn mass projection of Fig. 6 
due to /o(1500)r/ and /o(980)r?. In addition, there is some slowly varying contribution covering the whole Dalitz 
plots; it may come from the broad cr, i.e., /o(400 — 1200) in the Particle Data Tables [4]. We adjust fitted masses 
and widths of /2(1270) and 02(1320) by a few MeV from PDG values in order to achieve the optimum fits. This 
is because our main aim is to fit the production and decay angular distributions of these resonances. 






Figure 5: Data and fit (histogram) of invariant mass spectra for tt^tt^ (1 entry/event). 

Figs. 7 and 8 show differences on the Dalitz plot between fit and data. There arc small systematic discrepancies 
at the extreme right-hand edge of the Dalitz plot near an rjn mass of 1450 MeV. This discrepancy may be due to 
ao(1450) or 02 (1660) or /3(1405). The effect is small and cannot be analysed unambiguously into partial waves. 
Fits including these components have almost no effect on the main components of the fit, with the exception of 
rja, which covers the whole Dalitz plot and can absorb other small, ill-defined contributions. 

Figs. 9 and 10 show production angular distributions (after acceptance correction) for events lying in the 
/2(1270) mass band (1275 ± 100 MeV) and for events lying in the 03 (1320) mass band (1320 ± 50 MeV). It is 
immediately obvious that high orbital angular momenta are involved for both /2?7 and a27r at the higher beam 
momenta. The histograms show results of the partial wave fit described below. 



3 Formalism for Partial Wave Analysis 

For the tt^tt^?? final state, possible pp initial singlet states are 0~+, 2~+, 4 etc; for pp spin triplet, allowed 
states are 1"'""'", 2++, S"*""*", 4++, 5++ etc. For our case with center-of-mass energies below 2.41 GeV, only '", 
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Figure 7: Difference between Dalitz plots of fit and data where fit > data. 



2 '", 1^"*", 2++, and 4++ are expected to be significant [6] and this has been confirmed in our analysis; 
4 has been tried, but is not significant. The corresponding pp states with total angular momentum J, orbital 
angular momentum L and total spin angular momentum S in the usual contracted form ^^+^Lj are: ^5o for 
0-+, for 2-+, ■■'Pi for 1++, or for 2++, ^Fs for 3++, and or ^H^ for 4++. 

Let us choose the reaction rest frame with the z axis along the p beam direction. Then the squared modulus 
of the total transition amplitude is the following [16]: 
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where M is the spin projection on the z-axis in the initial state. The absence of M=0 for 1++ and 3+"'" is due to 
the vanishing of the Clebsch- Gordon (CG) coefficient (J = 2n + 1, Mj = 0|L = 2n+l, Ml = 0; S = 1, Ms = 0) 
with n as an integer. The relative minus sign for the interference term of {even)~^~^ and {odd)'^~^ partial waves 
with M=l and M=-l is also due to a property of CG coefficients. 

Each partial wave amplitude Ajpc includes contributions from various intermediate states (n), i.e., 



A-, 



CnAj 



(3) 



where C„ are free complex parameters to fit the data. In the present analysis, only /2(1270)?7, a2(1320)7r, 
ao(980)7r, (777, /o(980)77 and /o(1500)?7 intermediate states are considered. Amplitudes Ajpc^,^ arc constructed 
from relativistic Lorentz covariant tensors, Breit-Wigner functions and Blatt-Weisskopf barrier factors [17]. The 
amplitudes used for /o(1500)?7 and /2(1270)?7 intermediate states in our final fit are the following: 



G 



/o- 



(4) 



8 



- 6 





1 .05 GeV/c 



■■'■■tK':;; 





1 .525 GeV/c 



,ri;jjt :\ 



1 .642 GeV/c 




1 .94 GeV/c 



2 4 6 
m%nv) I GeV^/c' ) 



2 4 6 
Difference Dalltz plot (fit smaller than data) 

Figure 8: Difference between Dalitz plots of fit and data where fit < data. 
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where k^ is the four- momentum of the ry, G/„ = ~ ^-^tt ~ ''■^fo^ fo)~^ ^^'^ = i^h ~ ~ ^^h^S2)~^ 
are Breit-Wigner propagators for /q and /2. T^^ is a rank-2 tensor for /2 and is formed by the four-momentum 
(p) of /2 and its break-up four-momentum (q) as 
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Figure 9: Data and fit (histogram) of angular distribution da/d cosOr, for M^^r between 1175 and 1375 MeV (1 
entry /event). 
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Figure 10: Data and fit (histogram) of angular distribution da/d cosOtt for M^jr between 1270 and 1370 MeV (2 
entries/event). 
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The Blatt-Weisskopf barrier factors Bi{k) with a radius of 1 fm, the rank-Z tensors ig^^...s^ for pure i-wave orbital 
angular momentum of the r//o,2 system, and the spin- J wave functions (j)^^'"^^' (M) are standard as given in [17]. 

For /o(980)ry and arj intermediate states, the formulae are the same as for /o(1500)ry except for a different Gf^ 
for which we take the parameterization of Ref. [18], i.e., 

with Mr = 0.99 GeV, g„ = 0.117 GeV^, qk = 0.273 GeV^, m„ = 0.135 GeV and tuk = 0.496 GeV; 

Q = 1 + CqStttt ,201 

where Cq is a complex constant to be fitted by the data. 



r ( ] = r \/'^- '^'""'/' 1 + c.M\{s, - Ml)) 

"^'l + expiAiso-s))' ^1 - 16m2/M2 ^ ' 

with = 1.067 GeV, Gi = 1.378 GeV, P = 0.7 GeV, G2 = 0.0036 GeV, A = 3.5 GeV'^ and so = 2.8 GeV^. 

For aoTT and a27r intermediate states, the formulae are similar to those for forj and /2?7, but need symmetrization 
for two pions. The Breit-Wigner propagators for ao, a2, /o(1500) and /2 assume constant widths. The masses 
and widths (M, T) for ao and /o(1500) are fixed to be (0.9834, 0.085) GeV and (1.495, 120), respectively. Those 
for 02 and /2 are adjusted to fit the data. Based on these formulae, the data at each momentum are fitted by 
the maximum likelihood method. 

It is possible that the process pp r]w^w^ is driven, at least partially, by t-channel Regge exchanges. Even 
so, by Watson's theorem, each partial wave will acquire the phase variation of any s-ehannel resonance which is 
present; that is, amplitudes will contain singularities due to both s- and i-channel poles. Our strategy will be 
to express T matrices for individual partial waves Tl,j as sums over s-channel resonances. The formulae we use 
are 

ajPo^^is) = N^JAjPo^^{s)\\ (23) 

where s = M|p = M^^^, N is the normalization constant, ki and kn are the center-of-mass momenta of initial 
state and channel n respectively; Bl and Bi are barrier factors for the initial state and state n respectively; Anj 
are complex fitting parameters; Mnj and Tnj are masses and widths for resonances to be fitted. This prescription 
builds in the required threshold behaviour in each partial wave. By using a sum of resonances, we satisfy the 
constraint of analyticity. 



4 Results for partial waves 

The fit is shown as histograms in Figs. 5-6 for the mass spectra. It is obviously not perfect as regards broad, slowly 

varying components in the rnr projection of Fig. 6. However, since we are mainly interested in scanning the larger 
components from /2(1270)r?, arj, a27r and ao(980)7r intermediate states, we ignore those smaller contributions for 
the present study. 

The intensities of dominant partial waves are displayed in Fig. 11, and we shall discuss a fit to them below. 
The data points with error bars shown in Figs. 11 and 12 are our final fitted results for the partial wave cross 
sections ajpc^^ e&ch momentum for pp 7r°7r°?7 with r] — >■ 77. Small waves are displayed in Fig. 12. Partial 
waves with less significant contribution than those in Fig. 12 are dropped from our final fit. Table 2 shows the 
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Figure 11: Cross sections for partial waves making the largest contributions to pp — > tt^w^t] with ij 77. For 
diagrams with two components, the first label corresponds to the bigger component. The curves are the fit to 
the data points in the figure and the relative phases between components. 
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Figure 12: Cross sections for partial waves included in the final fit but giving smaller contributions than those 
in Fig. 11 to pp TT^TT^r] with rj — ?■ 77. 
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2++ 


2020 ± 50 


200 ± 70 


2.1 ±0.4 


4.3 ± 1.2 




2++ 


2240 ± 40 


170 ± 50 


2.5 ±0.6 


1.6 ± 1.6 




2++ 


2370 ± 50 


320 ± 50 


0.88 ± 0.64 


16 ±5 




1++ 


~ 1700 


~ 270 








1++ 


2340 ± 40 


340 ± 40 


0.6 ±0.6 


60 ±30 


0.84 ±0.53 


0-+ 


2140 ±30 


150 ± 30 


1.9 ± 1.7 


6±6 


10 ±5 


2-+ 


2040 ± 40 


190 ± 40 


3.0 ±0.3 


5.0 ±2.1 


0.4 ±0.2 


2-+ 


2300 ± 40 


270 ± 40 


2.8 ±0.7 


2.0 ±2.0 


0.5 ±0.5 



Table 2: Summary of fitted masses, widths and branching ratios corrected for their unseen decay modes. The 

mass and width of /4(2050) are fixed at PDG values, and the status of the 0~ state at 2140 MeV is questionable, 
as discussed in the text. The /i(1700) is beyond the accessible mass range. All states have I = 0, G = +1. 
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masses and widths of resonances included in the fit. Errors cover the range of values observed in a large variety 
of fits. The /i(1700) is below the range of masses accessible here, so its parameters are only approximate. 

The relative phases of the partial waves at each momentum are shown in Fig. 13. Since there is no interference 
between spin singlet and spin triplet, or between M=0 and M=l for spin triplet, there will be one overall phase 
undetermined for each M of spin triplet and for spin singlet. Hence we can only determine relative phases from 
our partial wave analysis. For spin singlet (0~ and 2~), the phases are relative to the partial wave of 2~ /2?7 
with L=0. For spin triplet with M=0, the phases are shown relative to 4+ — >■ a27r with L=3. For spin triplet 
with \M\ = 1, the phases are relative to 4+ — >• f2r] with L=3. 

4.1 = 4+ 

For 4++, a peak around 2090 MeV is clear for all 4++ channels. It can be fitted by a Brciit-Wigner amplitude 
with the mass and width fixed to the PDG values for the well established 4+ resonance /4(2050). The shift of 
the peak position to 2090 MeV is due to the centrifugal barrier factors for both initial and final states. Its decays 
into /2?? and a27r appear with comparable strength in the rjir^n^ channel. 

In addition to the /4(2050), there is clearly another 4++ peak around 2.32 GeV in 4+ — > /2r; in the M=l 
partial wave. This resonance may be identified with /4(2300) of the PDG, observed earlier in many analyses 
of pp — >■ 7r~7r+. The mass, width and phase with respect to /4(2050) are adjusted freely. The mass optimises 
at M = 2320 ± 30 MeV and the width at T = 220 ± 30 MeV. These agree closely with earlier determination 
quoted by the PDG, and also with recent VES data on ?77r+7r~ in the ttA reaction [19]. The latter find M = 
2330 ± lO(stai) ± 20{syst) MeV, F = 225 ± 20 ± 40 MeV. They also observe this resonance in data [20]. The 
/4(2300) is also observed in our data on pp — > tt'^tt" [15], with a slightly lower mass of 2295 MeV. The /4(2300) 
resonance acts as a valuable interferometer, determining the phases of S"*", 2+ and 1+ amplitudes over the mass 
range 2150-2400 MeV. 

From the M = 1 and M = amplitudes for 4+, we reconstruct the linear combinations for '^F^ and ^H4. 
Their intensities are shown in Fig. 14 for /2(1270)7? and a2(1320)7r channels. The /4(2050) resonance is almost 
purely ^i^4. The 027r channel is fed mostly by /4(2050) with a possible weak contribution from /4(2320); the 
contribution to a27r is barely significant. In contrast, the f2rj channel is fed by both /4(2050) and /4(2320) 
and the latter has a strong component. This is in agreement with the analysis of pp — >■ tt^tt"*" by Hasan 
and Bugg [6]; their Fig. 3 shows a strong component in /4(2320). The VES collaboration [19] finds that 
/4(2320) decays dominantly to /27?, in agreement with present results. 

4.2 = 3+ 

For J^*^ = 3++, there are significant enhancements at low mass (M ~ 2000 MeV) in both ao(980)7r and /2(1270)77 
with L = 1. At high mass {M ~ 2280 MeV) there is a strong peak in /2(1270)r; decays with both L = 1 and 
L = 3 decays. Fitted masses and widths are given in Table 2. There are no earlier listings of these resonances 
by the PDG. The observed phase with respect to /4(2050) and /4(2300) shown in Fig. 13 obviously requires the 
presence of at least one 3+ resonance, and is poorly fitted without two. The Argand diagram is shown in fig. 15. 

4.3 = 2+ 

For 2++, there is a peak in /2(1270)?? at ~ 2020 MeV and a peak at low masses in a2(1320)7r. At high mass 
around 2300 MeV, there is a strong peak in the a2(1320)7r channel. In /2(1270)?7, there is a further peak at 
~ 2230 MeV. The obvious question is how many resonances are required to fit these diverse structures. The 
phase variation observed on the Argand diagram. Fig. 15, requires at least two resonances from the observed 
360° phase advance. 

We find that the fit is poor without three resonances. The lowest peak fits naturally to a resonance with 
M = 2020 ± 50 MeV, F = 200 ± 70 MeV. Our data on pp -J> 7r°7r° independently find a resonance at 2020 MeV 
[15], and the analysis of Hasan and Bugg [6] of data on pp — )• 7r~7r+ likewise finds an /2 resonance at 1996 MeV. 
We have tried an alternative fit using instead /2(1920) observed by both GAMS [21] and VES [22] collaborations. 
The 2++ 0277 partial wave can be reproduced equally well with this assignment, but the 2++ — > f2r] partial 
wave is seriously underfitted by a factor 3 at 2050 MeV, ruling out a fit by /2(1920) only. 
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Figure 13: Relative phases (data points with error bars) obtained from the partial wave analysis and used for 
the fit (curves) to get Argand plots together with masses and widths of the resonances. The phases for 2+ and 
4+ with M=0 are relative to 4+ — )• a27r with L=3 and M=0; the phases for 1+, 3+ and 4+ with M=l are relative 
to 4+ f2'q with L=3 and M=l; the phases for 0~ and 2~ are relative to 2~ — >■ /2?7 with L=0. 
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Figure 14: Contributions to (a) /2(1270)77, (b) a2(1320)7r from (black squares) and ^H^ (open triangles). 

The proximity of this resonance to /4(2050) suggests that it may be identified as the qq state expected 
near this mass. Because 2+ amplitudes with M = 1 are negligible, and ^i^2 amplitudes have the same 
s-dependence; is the larger by a factor 1.44. This strong coupling of /2(2030) to pp also suggests 
identification with qq ^F2- high L in qq is likely to be associated with high L in decay channels, because of the 
peaking of wave functions at large r. 

At higher masses, a fit with a single resonance, shown by the dashed curve in Fig. 16, is much poorer than 
with two separate resonances. The peak at 2240 MeV in /2(1270)?7 has a mass compatible with ^(2230) observed 
in J/^ radiative decays [23], but has a larger width of about 170 MeV. This resonance may be interpreted as 
the n = 4 qq state, in the sequence /2(1270), /2(1565), /2(1920), /2(2240). The /2(2370) finds a natural 
explanation as n = 2 qq ^F2, i.e. the radial excitation of /2(2020). Its strong L = 3 decay supports this 
interpretation. In present data, both /2(2240) and /2(2370) appear in both ^P2 and ^F2, suggesting mixing 
between these states. 

The two peaks around 2020 MeV and 2370 MeV have masses and widths compatible with /2(2010) and 
/2(2340) listed by the PDG [4]. However, those observations were in the (jxp channel and could be different 
resonances, e.g. ss. We also remark that the peak in the 4>(j) data of Etkin et al. [24] actually appears at ~ 2150 
MeV. It is the rapid opening of the (jxf) phase space which leads to a pole at much lower mass, 2020 MeV, in the 
K-matrix fit to their data. 

4.4 = 1+ 

For J^^ = 1++ — > a2(1320)7r and ao(980)7r, there is a peak at the lowest masses. This suggests a resonance close 
to or below the pp threshold. However, as discussed below, the phase variation of the 1+ amplitude provides 
evidence for a resonance around 2340 MeV. The phase variation shown in Fig. 15 obviously requires resonant 
activity in the mass range 2000-2400 MeV. 

4.5 = 2- 

Partial waves with quantum numbers 2~+ and correspond to pp singlet states, and therefore there is no 
interference with even parity (triplet) partial waves. For 2 there is a strong peak in /2(1270)r? at ~ 2050 MeV 
and a smaller peak in a2(1320)7r at similar mass. There is evidence for a further peak at 2300 MeV. The 
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Figure 15: Argand plots corresponding to curves in Figs. 11 and 13. 
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Figure 16: Fits with /2(2240) (full curve) and without (dashed). 

lower peak is well fitted by a resonance with M = 2040 ± 40 MeV. The almost 360° phase advance observed on 
the Argand diagram points strongly towards the presence of two resonances, the second at 2300 ± 40 MeV. The 
PDG docs not list any "+ resonance in this mass range. A possible 1=1 partner is listed in the 

form of 7r2(2100). 

4.6 = 0- 

For there is a broad, slowly varying intensity with evidence for a strong peak superimposed at M ~ 2140 
McV. The slowly varying component may correspond to the broad object used in describing J/^* radiative 
decays to pp. lolu, K* K* ^ <j)(j) and rjirir [25]. The peak at 2140 MeV may correspond to a narrow resonance. 
However, it is observed in the rjcr channel, which contributes across the entire Dalitz plot. This contribution 
might absorb weak components not presently fitted to the data, for example due to ao(1450), 02(1660), p(1405) 
or further resonances in the production process around 2 GeV. In view of this possibility, the interpretation in 
terms of a resonance is ambiguous. Unfortimatcly. the relative phase with respect to 2~ is not well determined, 
so the phase variation cannot be used for independent evidence of resonant activity. 

5 Final fit to the partial waves 

To get more precise values for masses and widths for resonances, we use interfering sums of the Breit-Wigner 
amplitudes to fit the partial wave cross sections in Fig. 11 and the relative phases of the partial waves in Fig. 13 
simultaneously. The fit is shown in Figs. 11-13 as full curves. 

Besides the obvious resonances mentioned in the previous section, we need another 1++ resonance at about 
2340 McV with width ^ 340 MeV. Without it, we cannot describe the relative phase between 1++ and 4++ 
partial waves; also wc would need the lower 1++ resonance to be very narrow (< 50 MeV) in order to explain 
the sharply decreasing 1++ partial wave cross section. In our present fit with two 1++ resonances, the /i(2340) 
amplitiide interferes destructively with the tail of the lower 1++ resonance and causes the sharply decreasing 
cross section with a broad dip around 2340 MeV. The phase motion caused by this /i(2340) can be seen clearly 
in the Argand plots for 1++ partial waves of Fig. 15. 

In Table 2, the branching ratios are calculated at the resonance masses and are corrected for their unseen 
decay modes, except for ao(980) where Ta^^^ = Fog^^^^^. 
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For an ordinary qq state, the relative ratio /2fy/a27i' is expected to be smaller than 0.64. This allows for the 
36% component of ss in the rj. The centrifugal barrier and phase space will further suppress /2?7. Most of 
the branching ratios in Table 2 are in qualitative agreement with what is expected for qq states. However, the 

/2(2230) has an anomalously strong branching ratio to /2(1270)77 compared with a2(1320)7r. 

For the well-established /4(2050), only 44% of its branching ratios are listed in the Particle Data Tables [4], 

r - r 

in which nn has a branching ratio of (17 ± 1.5)%. In a very recent analysis [26] of pp tttt, the ratio "^i 

^ tot 

was reported to be (2.2 2.4) x 10^"^. Using this information, we can get the branching ratios of /4(2050) to 
pp, a2TT and f2V to be (1.4 ± 0.1)%, (30 ± 5)% and (3.9 ± 1.0)%, respectively. 

5.1 Comments on the resonance spectrum 

The /4(2044), /a (2000), /2(2020) and 772 (2040) cluster closely into a tower of resonances, as anticipated in the 
Veneziano model. Likewise the /4(2320), /a (2280), /2(2370), /i(2340) and r?2(2300) show indications of clustering 
into a tower at the higher masses. 

The 72(1920) originally discovered by both GAMS and VES has recently been confirmed in further VES data 
with increased statistics, decaying to ojuj [20]. There is also a strong /2 (1270)7? signal in VES 7y7r+7r~ data. 
Together with the /2(2020) we observe here, /2(2240) and /2(2370), this tentatively completes the identification 
of the qq I = and "^^2 states expected in this mass range. 

We conclude with some speculative suggestions of a scheme which concerns mixing of qq states with the 2+ 
glueball expected in this mass range. In our data on pp — )• 777777° [2], there is evidence for a further broad /2(1980) 
decaying to 7777, with mass M = 1980 ± 50 MeV, F = 500 ± 100 MeV. Its effects are seen clearly down to masses 
of ~ 1550 MoV. There is also evidence for a broad 2+ resonance in 4n final states in central production [27]. 
Such a broad state was predicted by Bugg and Zou [25]. It may be interpreted as a mixed state formed from 
the 2+ glueball, expected at ~ 2 — 2.2 GeV, and nearby qq states. Anisovich et al. [28] have argued that this 
mixing will lead to a broad state, accumulating the widths of nearby qq states and making them narrower. The 
/2(1920) and /2(2240) are indeed somwhat narrower this is usual for resonances in this mass range. Mixing with 
a glueball provides a natural explanation of the anomalous decays of /2(2020) and /2(2340) to (fxj), observed by 
Etkin et al. [24]. 

The glueball may be small, with radius ~ 0.3 fm; there are indications for this small radius in QCD Lattice 
calculations [29] . The small radius allows much of the glueball mass to be attributed to zero-point energy. Such 
a small object will mix preferentially with qq states rather than qq '^F2, whose wave; functions are strongly 
localised at large r. The preferential decays of /2(1920) and /2(2240) to /2(1270)??, despite its smaller phase 
space than a2(1320)7r, may be a further indication of mixing with the 2+ glueball. 

6 Summary 

In summary, we have observed a new decay mode rjirw for /4(2050). In addition, we have evidence for 7 new 
or poorly established resonances in the energy range from 1.96 to 2.41 GeV, i.e., /4(2320), /3(2000), /3(2280), 
/2(2240), /i(2340), 772(2040) and 7/2(2300). They appear to cluster into two towers of resonances around 2000- 
2050 MeV and 2300 MeV. Results are broadly consistent with earlier evidence for /4(2300), /2(2020) and /2(2340). 
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